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A natural generalization of the familiar second order Sturm-Liouville system is presented. This 
generalization consists of considering a number of differential equations defined on different intervals, 
instead of just one equation on one interval. The self-adjoint character of the differential equations is 
retained in the generalization, but the boundary conditions are relaxed considerably. The most general 
boundary conditions which can be accommodated by this sort of generalization of Sturm-Liouville 
theory are discussed. The existence of eigenvalues is proved, and a generalized orthogonality and a 
weak eigenfunction expansion theorem are derived. 
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1. Introduction 

The solution to a large class of problems of interest in applied mathematics can be reduced 
to the solution of the second order Sturm-Liouville system, 

(r(x)y' (x,k))'=(q(x)-\p(x))y(x,k), (1) 

where a ^ x ^ /;, and where 

ay(b, k) + 0/ (6, A) =yy(a, X) +8/ (a, X) =0. (2) 

The primes indicate differentiation with respect to x, and the proper choice of p(x) > 0, r(x) > 0, 
q(x), a, b, a, /3, y, and 8 is dictated by the particular problem. The theory of such Sturm-Liouville 
systems is well known [1, 2] 1 and will not be reviewed here. Suffice it to say that it is a very powerful 
technique for handling the problems to which it is applicable. 

Recently, generalizations of the second order Sturm-Liouville system have been developed 
which, while retaining the distinctive Sturm-Liouville character, are applicable to problems 
involving composite media [3, 4, 5, 6]. The purpose of this paper is to discuss a different general- 
ization of the second order Sturm-Liouville system which retains the Sturm-Liouville character. 
In particular, an eigenvalue equation, an orthogonality theorem, and a weak eigenfunction expan- 
sion theorem will be derived. This eigenfunction expansion theorem is the first place where the 
results of this paper deviate from the previous work on Sturm-Liouville systems. It will be shown by 
counter-example, that the eigenfunctions of this Sturm-Liouville like system do not form a com- 
plete set in the sense of the previously described Sturm-Liouville systems. 

The generalization of the second order Sturm-Liouville system discussed in this paper, and 
that which is applicable to composite media problems are similar in many ways. Each arises when 
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attempting to solve time dependent heat transfer problems in complicated systems, and each 
consists of considering a number of differential equations defined on different intervals, instead 
of just one equation defined on one interval. In both generalizations the self-adjoint character of 
the differential equations is retained, but the conditions which the solutions must satisfy at the 
ends of the intervals are relaxed. In both cases, the solutions to all of the equations are coupled 
by relations involving their values and the values of their first derivatives at the ends of the intervals 
on which they are defined. The difference between the two generalizations is in the specific form of 
the coupling. In the case of the composite media generalization, the coupling connects the boundary 
values at the right end of the first interval with those at the left end of the second interval; the 
boundary values at the right end of the second interval with those at the left end of the third interval, 
and so forth. In the alternate generalization, introduced here, the coupling connects the boundary 
values at the left ends of all of the intervals. Figures la and lb show examples of the two types of 
coupling schematically. It is interesting that the form of coupling introduced here results in a 
simpler eigenvalue equation than that which results from the composite media generalization. In 
the case where only two intervals are considered, the composite media generalization is a special 
case of the generalization introduced here. Except for this case, this generalization of the second 
order Sturm-Liouville system does not appear to have been studied previously. An example of an 
application of the theory developed in this paper, is the calculation of the temperature distribution 
as a function of time in a number of wires whose left ends are all welded together and whose right 
ends are all connected to a heat sink. 



B 





FIGURE 1. Schematic representation of different types of Bound- 
ary value coupling 

— Interval on which self-adjoint differential equation is defined. 

° End point of interval where Sturm-Liouville boundary conditions are imposed. 

• End point of interval where boundary conditions appropriate to the composite media 

generalization are imposed. 
O End point of interval where boundary conditions appropriate to the generalization 

introduced in this paper are imposed. 



2. Postulates and Principal Results 

Consider the N self-adjoint differential equations obtained by letting i= 1, 2, . . ., N in 
the equation, 

(n (xi )y'i (xi,k))'=(qi (xi ) - kpi (xi ) ) y t (xi , A ) (3) 
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which is to hold for a* ^ x% ^ bi, where r%(xi ) > and pt(xi ) > 0. For the present, retain the usual 
Sturm-Liouville boundary condition, 

aw(h,X)+j8<ri (fei,X)=0, (4) 

at the right end of each interval, but at the left end of the intervals employ a generalized bound- 
ary condition, consisting of the N—l equations obtained by letting i=l, 2, . . ., TV— 1 in the 
equation, 

yiyiioiX) + hiyl (a>i^) = yi+iyi+\ (a i+u k) + S i+1 y^ x (a i + 1 , X) ^ 0. (5) 

One more coupling condition is needed to give a total of 2N boundary conditions, so introduce 
the general relation, 

y 2Ciy i (ai,k)+D i y , i (a i ,\)=0, (6) 



Ay,-C/o,>0, (7) 



for 1=1,2, . . . , N. 



Equations 3, 4, 5, 6 and 7 define the system. Now it is necessary to determine under what con- 
ditions a set of functions, (y, (%i,X), y-> (#2, X), . . . , y.\(x\,k)), exists such that these five equa- 
tions are satisfied. 

First concentrate attention on eqs (3) and (4). For any positive integer value of i ^ N, it can 
be shown, following Ince [7], that for every value of X, eqs (3) and (4) determine, except for the 
specification of a multiplicative constant, one and only one solution which is continuous and 
has a continuous first derivative on the interval, (a* ,6*). Keeping this in mind, divide eq (6) by 
Xiyi (a\ , X) + S\y' {a x , X) and substitute the TV — 1 equations represented by eq (5) into the result, 
to obtain the equation, 

" Cm (a h k)+D iy ; (a,, X) 
£ ytyi (auk) + 8 iy ! (a u k) U ' w 

which is generally a transcendental equation in X, and which is independent of the multiplicative 
constant which is as yet unspecified in y- x (xu X). Equation (8) is called the characteristic or eigen- 
value equation of the system, since its roots are the eigenvalues of the problem; that is, the only 
values of X for which eqs (3), (4), (5), and (6) are mutually consistent. The existence of a countably 
infinite number of such roots is proved in appendix A of this paper. Let k m be the mth root of eq 
(8), then it is called the mth eigenvalue of the system, and yixuk m ) will be called the mth eigen- 
function on the ith interval of the system. This completes the discussion of the eigenvalue equation. 
Now attention will be turned to exhibiting the Sturm-Liouville like properties of the eigen- 
functions. Other than the mere existence of eigenvalues, the most distinctive characteristic of 
Sturm-Liouville theory is the theorem expressing the orthogonality of the different eigenfunctions 
and the related eigenfunction expansion theorem. For the system described in this paper, the 
orthogonality theorem, which is proved in appendix B, takes the form, 

N rbi 

2 Wi Pi(xi)yi(xu k m )yi(xi, k H )dxi = 0, (9) 
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for m # n, where 



D,y,-C& 
Wl n(a») (10) 



Now the problem of an eigenfunction expansion theorem will be considered. Since the func- 
tions, y%{xu ^m), are determined only to within a multiplicative constant, it is convenient to define 



dim 

and 



fbi 11/2 

Pi(xi)yi(xi, \ m )yi(xu km)dXi\ , (11) 



<t>im(Xi) = yi(xu k m ) I dim. (12) 



Thus each of the functions, <f)i m (xi), is square normalized to unity on the ith interval with respect 
to the weighting functions, Pi(xi), and this property uniquely determines the function, <f>i m (xi). 
Of course, the dim are as yet undetermined, since the yi(xi, k m ) are determined only to within a 
multiplicative constant. 

Now consider the possibility of choosing these multiplicative constants (and hence the dim) 
in such a way as to express the A^ arbitrary functions, gi(xi), 1=1, 2, . . ., iV, as a sum of 
eigenfunctions. For any given value of i, the constants, dim, which are obtained by varying m are 
independent of each other. However, when m is held constant and i is varied, the resulting dim 
are not independent. Once one, say dj m , is specified all of the others can be obtained from the 
equation, 

d im ^ 7j <fon (aj) + 8j <f>jm ( Oj ) (J3) 

djm 7i<l>im(ai) + SitfimiOiY 

which is obtained by substituting eq (12) into eq (5). Thus in the eigenfunction expansions, 

00 
fi( x d = ^ dim<j)im(Xi) , (14) 

m=\ 

the TV functions, fi(xt), cannot be chosen arbitrarily. Once one of them is specified, all of the rest 
of them are uniquely determined. This is a strong restriction, but it does not necessarily prevent 
the theory from being useful. However, it does weaken the eigenvalue expansion theorem con- 
siderably. The theorem stated below and proved in appendix C is the strongest that can be proved 
when the functions to be expanded are left completely arbitrary. 
For each value of; = 1, 2, . . ,,N, the function, 

X oc 

M X J) =2 d Jm<l>jm(Xj)=^ yj(Xj, km), (15) 

»t=l m=l 

will be a "best" approximation to the function, gj{xj), in the sense that the expression, 

I=t Wi I"' Pi(xi)[gi(xi)-Mxi)rdx u (16) 

i=i Ja i 

is a minimum, if and only if for each m= 1, 2, . . ., 

N Cb { 

^RmjiWi I Pi(Xi)gi(Xi)(t)im(Xi)dXi 

djm=- J (17) 

J,(R mji )Wi 

i=\ 
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where 

R _ dim Jj^jm ( dj ) + dj^pn ( CLj ) 

o/wi y/0/;» ( a,- ) + Oi0,- m ( a, ) 

Since the above eigenfunction expansion theorem does not require that each function, ft(xi) , be 
identical to the function, gi(xi) , a criterion for goodness of fit is useful. It is shown in appendix 
C that the ratio, 

R = 2 2 *w») 2 / £ Wi f w(«i) fe(^)] 2 ^ d9) 

n=i ?=i ' i=i ^ "/ 

has the property that, ^ R ^ 1, and R = 1 if and only if /i(#i) = gi(#i), for i= 1, 2, . . ., N. Thus 
/? can be used as a criterion for goodness of fit in a qualitative manner. Furthermore the set of 
functions for which the eigenfunctions of the problem form a complete set are characterized by 
R=l. 

3. An Example 

As an example, let £ = 1, 2, and 3, in eqs (3), (4), (5), and (6), and let a\ = — 1, bi = 0, e=X 1/2 , 
ri(xi) = l, qi(xi)=0, pi(xi) = (i) 2 , a*=l, j8i = 0, yi=l, 8f = 0, Ci = 0, and D< = 1. The resulting 
system arises when attempting to determine the temperature distribution, Ti(xi, t) , as a function 
of time in three straight wires of identical thermal conductivity and area, whose specific heats per 
unit volume are in the ratio 1, 4, 9, and whose right ends are attached to a block whose thermal con- 
ductivity and heat capacity are infinite, and whose left ends are attached to a block whose thermal 
conductivity is infinite, but whose heat capacity is zero. If the initial condition for this problem is 
that Ti(x- t , 0)=— xu for i= 1, 2, 3, then T\{x\, t) can be written in the form 



Ti(x u = 2 <hm<l>im(xi) exp [ - (e m ) 2 t] , (20) 



where €m=X w is the mth root of eq (8), if and only if fi(xi) = — #i,for i= 1, 2, 3, where fi(xi) is 
defined in eq (15). However, this is not the case. The functions, 



Zi(xi) = ^yi{xu\ m ), (21) 



are shown in figure 2. It is clear from the figure, that these functions are not very good approxima- 
tions to — x\. For this particular case, R = 0.990159. There is no significant difference between 
these results and the results when only the first fifteen terms are included in the sum, and presum- 
ably there is no significant difference between these results and the results when all of the terms 
are included in the sum. Thus in order to obtain Tt(xu t) it is necessary to add to the right hand 
side of eq (20), a nonseparable function of x, and t, which is a solution to the corresponding time 
dependent heat transfer equation for the problem. The properties of such non-separable solutions 
remain to be investigated. 

Now suppose that we eliminate the third interval in the above problem, restricting i to be 1 
or 2. Then the problem is equivalent to a composite media problem, and interestingly enough the 
eigenfunctions form a complete set as can be seen from figure 3. In this case, /? = 0.999989, for 
just the first fifteen terms of the sum. 
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Figure 2. Graphs o/Z, (x ( ) for i= 1, 2, and 3, showing that Z, (x ( ) is not a very good approximation 
to —xi in the three interval example considered in the paper. 



4. Discussion 

The fact that nonseparable solutions appear in a prominent way in the generalization of Sturm- 
Liouville theory discussed here certainly limits the usefulness of the theory to some extent. How 
serious this limitation is in various applications will depend upon the general properties of these 
non-separable solutions. Thus an investigation of such non-separable functions would appear to be 
quite worthwhile. 

Nowhere in the derivation of the results presented in this paper does the assumption that 
a; < xi < bi play a fundamental role. All of the results would have been identical if for some of the 
values of i, bi < xi < a,-, but the discussion would have been somewhat more complicated. To further 
generalize the second order Sturm-Liouville system, the boundary conditions at the points, bi, 
could be relaxed. It is possible to replace them by generalized boundary conditions identical in 
form to those stated earlier in this paper, by introducing a number of new intervals on which more 
self-adjoint differential equations are defined. The type of coupling envisioned here should be 
obvious after glancing at figure lc. The detailed development of this further generalization of the 
second order Sturm-Liouville system introduces no new ideas, and it is fraught with notational 
difficulties. Consequently it will not be pursued here. However, a few remarks are in order. 



126 



CD 



1 = 1 
1=2 




"iToo 



Figure 3. Graphs o/Z,(x,) for i = l, and 2, showing that Z,(x,) is an excellent approximation to 
— Xj in the two interval {composite media) example considered in the paper. 

First, cross-linking among the branching intervals will not generally allow a nontrivial solution. 
But it is possible to cross-link any two intervals for which r(a) = r(b) for both intervals. Second, 
both the composite media generalization and the generalization described earlier in this paper are 
special cases of this generalization of the second order Sturm-Liouville system. Third, in a detailed 
development of this theory, it will not be possible to derive a simple eigenvalue equation as was 
done here. A determinantal equation, similar in principle to, but more complicated in actual detail 
than that in reference [4] will have to be derived. 

5. Appendix A 

In this appendix, eq (8) will be proved to have an infinite number of roots. In discussing eq 
(8), it is convenient to define, 



Fi(k) 



and, 



Ciyi{ai,k) + Diy'i («,-,\) 
> ( y l (a f ,\)+8,y, ! (ai, X) ' 



F(X) = TF«(X). 



(1A) 



(2A) 
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First it will be proved that for each i = 1 to N, F;(X) has an infinite number of singularities, and 
that it is positive in some interval immediately to the right of each singularity, and negative in 
some interval immediately to the left of each singularity. In order to prove this, it is necessary 
to consider each of the four cases, C/ = 6\ = 0, Z); = y; = 0,C;yj =N0, and DiSi 4= 0, separately. Every 
possible choice of C/, y>, Du and 8i which is consistent with eq (7) is covered by at least one of 
these cases. Before looking into each of these cases separately, it is necessary to review some 
of the properties of yi(ai,\). Ince [8] has shown that if yt(xi,\) satisfies eqs (3) and (4), then the 
equation, yi(ai,k) = 0, has a countably infinite number of roots. These will be denoted by /uiim, 
where the /Jiim are labeled so that film < (J<im + i. Ince [9] has shown that as X is increased from 
(JLim to /Xi m +i, the expression, 

yi{(H,k) 

decreases monotonically from + 00 to — °°. Thus y f t (ai,k) also has an infinite number of roots, 
denoted by co* m , and as X is increased from cotm to a>im + i, the expression, 

y[{ai,\) ' 
increases monotonically from — 0° to + °° . 

Case 1: Ci = 8i = 
In this case, 

Diy'i(ai,\) 



F i {\) = 



yin(a u k) ' (3A) 



and dividing eq (7) by (y0 2 shows that A/y, is a positive number. Thus Ft(k) decreases mono- 
tonically from + 00 to — °°, as X is increased from film to /Ltim+i, so it has an infinite number of 
singularities and it is positive immediately to the right of each singularity and it is negative im- 
mediately to the left of each singularity. 

Case 2: A=7i = 

In this case, 

f;(x)= £y^Ai (4A) 

but dividing eq (7) by (S;) 2 shows that Ci/8i is a negative number. Thus Fi(k) decreases monotoni- 
cally from -foo to — °° as X is increased from coim to Q)tm + i. So again, F%{\) has an infinite num- 
ber of singularities, is positive to the right of each singularity, and is negative to the left of each 
singularity. 

Case 3: 0(8,4=0 
\Di yi{a u \)J \8i yia u k) ) 



then F 'i(X) =(Dil 8 i)fi(k). Now the equations, 



y!(ai,X) Ci 
yi(auk) D t 
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and 

y-(gjA) y, 
yi(auk) 8j 



= 0, (7A) 



each have one and only one root between i±\ m and fJLim+1- This root of eq (7A) will be denoted by £, m 
and the equivalent root of eq (6A) will be denoted by i)i m . First supposing that < Dfiu use eq (7) 
to show that — -y,78/ < —CJDu so that jjL im < r)i m < tjim < fiim+i. Thus 

— Yi yfauk) C, 
8i ywuk) D, (8 A) 



and 



y'iiauk) ji 

— TT < — T for f im < A < (A 



7 ( — \\ < ~~R ior ^ hn< X< i^'w+i- (9 A) 

Since £, m is a root of eq (7 A), //(A) has a singularity at £,- m , and eq (8 A) can be used to show that //(A) 
is negative when 17 | n < A < f,- m , and eq (9A) can be used to show that //(A.) is positive when £,-,„ 
< k < fJLim+i* Letting m= 1,2, . . . in the above discussion shows that //(A) has an infinite number 
of singularities and that it is positive immediately to the right of each singularity and that it is neg- 
ative immediately to the left of each singularity. Now, supposing that D,8, < 0, use eq (7) to show 
that 



' C «.<-2! 



Di 8, 



SO that flim < gim < Tfim < fJLim+l. 
Thus 



d y/(a,-,A) . yi 
-7r<— — -<-— forf/, w <A<T7,- m , (10A) 



and 



y> y/(a/,A) 
~8" i < ^(^X)' for/A "" <X< ^» ! ( 11A ) 

Since gi m is a root of eq (7A),/) (A) has a singularity at £,- m , and eq (10A) implies that// (A) is negative 
when £,,„ < A < 77^, and eq (11 A) implies that /J (A) is positive when /x im < A < £i, ft . Letting m= 1, 
2, ... in the above discussion shows that //(A) has an infinite number of singularities and it is 
negative in some interval immediately to the right of each singularity and that is positive in some 
interval immediately to the left of each singularity. So in both the case where Did-, > and in the 
case where A8? < 0, Fi(k) = (Dt/Si) fi(k) has an infinite number of singularities and it is positive 
immediately to the right of each singularity and negative immediately to the left of each singularity. 

Case 4: C iyi 4= 

The proof is similar at each step to that of Case 3 and will not be presented here. 
If Fi(k) has a singularity at K,- m , then eq (2A) implies that F(k) also has a singularity there, 
and F(k) behaves like F,(A) in some neighborhood of Kim* Thus, if Km and Kj n are two consecutive 
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singularities of F(X), and Km < Kj n , then F(X) is positive immediately to the right of Kim and 
negative immediately to the left of k>. Now since F(\) is continuous for Kim < X< Kj n , it must 
have a zero in this interval. Thus eq (8) has at least one root between any two consecutive singu- 
larities of F(\). Finally, since F(X) has an infinite number of singularities, it has an infinite number 
of roots, as was to be proved. 



6. Appendix B 

In this appendix, eq (9) will be proved. The proof requires a lemma which will be proved 
separately, but in somewhat more general form than needed. 

Lemma 

Since this lemma consists almost entirely of algebraic manipulations, it will prove convenient 
to use the shorthand notation, 

Mi = Mi(ai,£), (IB) 

Ni=Ni(a u r)), (2B) 

and 

ri = n(ai), (3B) 

where M;(x,, £) and Ni(xi, 77) satisfy eqs (3), (4), and (5), with X= f , and k = 7), respectively. Now let 

S = ( yi Ni + SiN't ) 2 (CjMj + DjM'j ) - (ytMi + S,-M; ) f (CjNj + DjN' } ) , (4B) 

and notice that if Mj and Nj both satisfy eq (6), then S = 0, but if only one of them satisfies eq (6), 
then S * 0. 

Since both Mj and Ni satisfy eq (5), 

S=f [{CM^D^iy^^b^-iC^+D^iyMs^m^ ' (5B) 

j=i 

When the indicated multiplications in eq (5B) are carried out, some of the terms cancel one another, 
and the remaining terms can be combined to yield, 

S= f Wm-CA) (M^Nj-N^Mj) • (6B) 

Thus by eq (10), 

S^fWjrjWNj-N'M) (7B) 

Churchill [10] has shown that eqs (3) and (4) are sufficient conditions for 

(€-v) r i pAxj)Mj(xj,{)Nj(x j ,r ) )dx j = r j (MM-N;M j ). (8B) 

Jaj 
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Thus 

(f—q) f Wj \ bj Pj(xj)M j (x j ^)N j (xj,7 i )dxj = S 9 (9B) 

where S is given by eq (4B). This is the significant result of the lemma. 

Now let €=k m and 7) = \ n - Then Mj(xj, £)=yj(xj, k m ) and Nj(xj, T)) = yj(xj, X w ), so Afj and 
Nj both satisfy eq (6), so S = 0. And if m//i, eq (9B) becomes eq (9). This completes the proof. 

7. Appendix C 

In this appendix, it will be shown that eq (17) is a necessary and sufficient condition for the 
expression, 7, defined in eq (16), to have a minimum with respect to variation of the djm subject to 
the constrants imposed by the N — 1 equations represented by eq (13). 

First notice that by the restrictions stated after eq (3), by eq (7) and eq (10), 7 ^ 0, so it has a 
minimum. Also notice that 7=0, if and only if fi(xi) — gi{xi) for £=1, 2, . . ., N. Next take the 
variation of 7 with respect to djm where m=l, 2, . . . , and j is fixed. Using eq (13) to justify re- 
placing the partial derivative of dim with respect to djm by di m /dj m , rearrange the result to obtain, 

S(/)= f 8(d jm ) Y ($*) Wi I"* Pi (x i )( gi (x i )-f i (x i ))ct> i „ l (xi)dx i = 0. (1C) 

Now since each dj m is independent of dj n for n 4= m, the coefficient of every d(dj, t >) must be zero. 
Combining eqs (15), and (18) with this result yields 

n n n i N Cb { 

V RmjiWi Pi(xi)gi(xi)<t> im (xi)dxi = -7-X Wi Pi(xj) fi(xi)yi(xi, \ m )dxj, (2C) 

for 771=1, 2, . . . . Finally applying eqs (15), (9), and (11) to the right-hand side of eq (2C), and 
again using eq (18) one sees that the simultaneous equations are diagonalized and that they are 
equivalent to eq (17). All of the steps in this proof are reversible, so eq (17) is a necessary and 
sufficient condition for I to have a minimum value as was to be shown. 

In order to obtain a criterion for goodness of fit, expand I to obtain, 

/=fr, f* Pt(*t)<a(*))*dxt+f, Wi f%,(*i)(/ite)) 2 <fe-2 Y Wx I"* Pi(xi) gl (xi)fi{xi)dxi. 

1=1 J"i (=1 J"i i=l J"i 

(3C) 

Notice that if eq (2C) is multiplied by dj m and summed over m, the result is 

2 W t \ bi Pi(xi)gi(xi)Mxi)dxi= V Wi I' 1 pt(xi){fi(xi))*dxi. (4C) 

i=i J*t i=i J a i 

Substitute this result into eq (3C) to obtain 

/= S W * r Pi(xi)(giM) 2 dxi-y,Wi { bi Pi(xi)(fi(xi))*dxi. (5C) 

Finally substitute eq (15) into eq (5C), and use the orthogonality theorem, eq (9) to obtain, 

1=2 Wi I" 1 P i(xi) ( gi (xi)Ydxi- 2 f Wi(d im )\ (6C) 

i=i J a i m=i i=i 
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Since 7^=0, 0=^/?=^l, where R was defined in eq (19). The second property of/? follows from the 
property of I that it is zero if and only if gi(xt) =fi(xi) for a* ^ x% ^ bu and i= 1, 2, . . ., N. 



I thank R. C. Kraft and W. P. Reid for constructive criticism on various aspects of this work, 
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